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Abstract 



Most models of inflation have small parameters, either to guarantee suffi- 
cient inflation or the correct magnitude of the density perturbations. In this 
paper we show that, in supersymmetric theories with weak scale supersym- 
metry breaking, one can construct viable inflationary models in which the 
requisite parameters appear naturally in the form of the ratio of mass scales 
that are already present in the theory. Successful inflationary models can be 
constructed from the flat-direction fields of a renormalizable supersymmetric 
potential, and such models can be realized even in the context of a simple 
GUT extension of the MSSM. We evade naive "naturalness" arguments by 
allowing for more than one field to be relevant to inflation, as in "hybrid 
inflation" models, and we argue that this is the most natural possibility if in- 
flaton fields are to be associated with flat direction fields of a supersymmetric 
theory. Such models predict a very low Hubble constant during inflation, of 
order 10 3 -10 4 GeV, a scalar density perturbation index n which is very close 
to or greater than unity, and negligible tensor perturbations. In addition, 
these models lead to a large spike in the density perturbation spectrum at 
short wavelengths. 
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I. INTRODUCTION 



Inflationary models M in general require small parameters in the particle theory La- 
grangian, to provide the flat potential needed for sufficient inflation and for the correct 
magnitude of density fluctuations. The need for unmotivated small parameters tends to 
weaken the credibility of a theory, so one hopes that the origin of these parameters can be 
understood. It is conceivable, of course, that the explanation lies beyond our present under- 
standing, just as we presently have no accepted explanation of why the Yukawa coupling of 
the electron is 2 x 10~ 6 , or why the weak scale lies 17 orders of magnitude below the Planck 
scale. Nonetheless, it would be encouraging to find that the small parameters required by 
inflation could be obtained from small parameters that are already essential to the particle 
theory, so that no additional small parameters are introduced. Models in which the small 
parameters arise as ratios of known particle physics mass scales are particularly attractive 

From a field theoretical perspective, it is difficult to see how flat direction fields can be 
present in a nonsupersymmetric theory (with the exception of Goldstone bosons, considered 
in Ref. ||) given the effect of radiative corrections. We will assume therefore that the 
world is super symmetric, and ask whether an inflationary potential can arise naturally in 
the context of the mass scales which we expect might be present. Some natural candidates 
for these scales could be the Planck scale M p « 10 19 GeV, the GUT scale Mq « 10 16 GeV, 
the intermediate scale Mj 10 11 GeV, and the supersymmetry breaking scale 7713/2 ~ 1 
TeV. 

In the context of supersymmetric models, an attractive scale for the vacuum energy 
density during inflation would be set by the intermediate scale, Mj. This scale is very likely 
to be present in a hidden sector model of supersymmetry breaking. It is also the right energy 
scale for the potential associated with moduli fields, which might be natural candidates for 
flat directions. The problem here, however, is that simple dimensional analysis arguments 
(to be reviewed in Sec. II) show that density perturbations would generically be either 
far too small 0|| or far too large, depending on assumptions. For the case of inflation 
driven by a single moduli field, the dimensional arguments show that the requirements 
of sufficient inflation and correct density perturbations imply that 1) the variation of the 
inflaton (moduli) field during inflation is of order M p , and 2) the energy density during 
inflation is of order M 4 where M ~ 10 16 GeV. An energy density of order Mf would 
produce density perturbations too small by about ten orders of magnitude. For the case of 
a chaotic inflationary scenario, the variation of the inflaton is again of order M p , but in this 
case the density perturbations are much too large unless the quartic coupling A is about 

10- 12 . 

In this paper, we show that the argument that inflation at an intermediate scale is 
untenable lacks sufficient generality, and can evaporate if one drops the assumption that 
inflation is driven by a single scalar field. We describe a class of two-field models, for which 
dimensional analysis estimates show that 1) the variation of the inflaton is of order Mj or less, 
and 2) the energy density during inflation is of order Mf. We then go on to illustrate these 
ideas with models motivated by supersymmetry with soft supersymmetry breaking. We will 
find that these models not only solve the naturalness problem of obtaining sufficiently many 
e-foldings of inflation, but also generate very nearly the correct size of density perturbations 
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based on the parameters of supersymmetry breaking. We therefore refer to our models as 
"supernatural" inflation. 

This model contains a similar structure to the "hybrid" inflation models, proposed by 
Linde and studied by Copeland, Liddle, Lyth, Stewart, and Wands [|]. The fact that the 
standard dimensional naturalness arguments for the number of e-foldings and for Sp/p do 
not apply, and that the Hubble scale during inflation will be low was also clearly recognized 
by these authors. Our point here is to emphasize that the most natural scales for successful 
implementation of two field inflation of the "waterfall" type are the scales associated with 
supersymmetry breaking and the Planck scale. Furthermore, our models more accurately 
reflect masses and couplings associated with flat direction fields, and we will motivate the 
parameters and potential we use by consideration of flat directions in the MSSM. Hybrid 
inflation in the context of SUSY leads one to the interesting conclusion that the Hubble 
scale during the inflation which established the density perturbations might have been of 
order 10 3 -10 4 GeV, rather than 10 13 GeV. 

In the following section, we present the general arguments for why supersymmetry scales 
do not work in single field inflation models. We then review the general idea of "hybrid" or 
"waterfall" || models, and show why the single-field arguments do not apply to the two-field 
case. In Sec. Ill, we present supernatural inflation models, in which we assume the inflation 
sector consists of flat direction fields whose potential is generated through supersymmetry 
breaking and operators in the superpotential which do not permit the two fields to be 
simultaneously flat. We derive the constraints on parameters consistent with the requisite 
number of e-foldings and density perturbations. We conclude in the final section. 

II. ONE VS. TWO FIELD INFLATION 

We begin this section by reviewing the "standard" arguments for why the inflaton in 
"natural" inflationary models varies on the scale M p and why the scale for the energy density 
should be larger than the intermediate scale in inflationary models with a single field. 

For the purposes of these dimensional arguments, we first assume the potential takes the 
form 

V = M 4 £(0//) , (1) 

where Q is a bounded function of order unity. Here we have in mind for example a moduli 
field, with M « Mj. If we assume the slow roll equation of motion 3H<p = —V, where H is 
the Hubble constant during inflation, the number of e-foldings is 

J J V <p J V M%V M p M p Q' v ; 

There are essentially two possibilities. If Q is a bounded function, and Q' is not very tuned 
to have very flat sections, one is in the regime of what might be expected for a moduli type 
field. In this case, the requirement of about 60 e-foldings of inflation favors / of order M p 
and a change in <fr during inflation at least of order M p . Even when this is satisfied, some 
tuning of the potential is required. 
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The alternative possibility is that one is in a chaotic inflationary scenario, in which 
case Q will be dominated by monomial behavior for sufficiently large field, and V ~ V/<f>. 
In this case, / is not defined, but one would still conclude A0 rs M p . 

Density fluctuations are also readily estimated under the assumed form of the potential. 
They are given by 

P ~ ~ V ~ \M p ) M p Q> ■ U 

Assuming a potential of the moduli type, with Q and Q 1 of order unity and / of order M p , 
we find that Sp/p ~ (M/M p ) 2 favoring M ~ 10~ 3 M P . Detailed calculations might change 
M by an order of magnitude or so, but it is clear that M ~ 10 11 GeV ~ Mj is strongly 
disfavored. 

In a chaotic scenario on the other hand, one would conclude that the density fluctuations 
are too large unless there is a small parameter. For example, a simple dimensional argument 
would lead to the conclusion that for V = A0 4 , A ~ 10~ 12 . Without further motivation for 
these small numbers, such a potential seems unlikely. 

So one is led to the conclusion that it is difficult to naturally obtain sufficiently many 
e-foldings and the correct magnitude of density perturbations, without invoking either small 
numbers or a new mass scale. 

It is apparent, however, that there is a loophole in the above argument. From Eq. (^) 
it is clear that the constraints on / and A0 during inflation arise because it is the same 
potential V(0) that controls the inflation rate H and the speed of the inflaton field 0. These 
constraints can be avoided, therefore, if the energy density during inflation is provided from 
some source other than the scalar field which rolls and controls the ending of inflation. 

The simplest way to implement this idea would be with two fields. This idea is essentially 
that first proposed by Linde |J as "hybrid" inflation or "waterfall" models. There are two 
fields if> and 0. The first field, which we call the inflaton, has a very flat potential. It starts 
at a large field value, and slowly rolls (via its classical field equations) to the origin. 

The second field, 0, has a potential whose minimum is far from the origin. In most 
previous incarnations of hybrid inflation, the scale of variation of this field is Mj, though 
in our models the scale will be M p . When if) has large field value, it gives a positive mass 
squared term in the potential at the origin, so the classical field equations push to the 
origin. When if) gets sufficiently small (of order Mj or less in our models), the mass squared 
of goes negative, and makes the transition from the origin to M p . 

The key feature of this model is that the energy density during inflation is dominated by 
the potential energy of the field at = 0. There are no tunings in the if) potential to get 
a small mass during inflation and a large mass afterwards, since its mass is always small, as 
is its potential energy. Because H depends on the value of and is not determined by the 
field if), which acts as a switch to end inflation, the naive estimates do not apply. 

The second key feature of this model is that the ending of inflation is controlled by when 
the mass squared at the origin changes sign. One can obtain a large number of e-foldings 
with the variation of the inflaton field if) much less than M p . Let us see this explicitly. 

We assume a potential which takes the form 

V = M 4 £(|0|//)+s(|0|,H) + m 2 H 2 , (4) 
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where M and / are to be determined, the function g is the term responsible for the if> 
dependence of the (f> mass, and m is of order 
We now have 

N „I^„-J^*J£«). (5) 
J ip J m 2 ip m 2 M 2 yPfinal) 

Notice that the scale M in the numerator is independent of the mass and coupling of the ip 
field (in the limit that the ip contribution to the energy density is small) so that the previous 
arguments for one-field inflation no longer apply Clearly for inflation to give several e- 
foldings requires only that ip changes by an order of magnitude, and that M 4 > m 2 M 2 . No 
inflaton variation of order M p is required, and so far, it seems M ~ Mj could be a good 
choice. 

Let us now consider density fluctuations under the same assumed form for the potential. 
We find 

5p H 2 H 3 M 6 

p tp m 2 ip M^m 2 ip 

The point is that the numerator H 3 has its scale set by the <p potential energy while the 
denominator is determined by the ip field. We construct a model so that ip a t the end of 
inflation is of the order Mj or smaller. If we also take M ~ Mi, we find Sp/p ~ Mj/M p 
or bigger (rather than (Mj/M p ) 2 as was the case in single field models). Although the 
coupling between the fields can have a coefficient which varies by many orders of magnitude, 
as does ip in Eq. (||), the strong M dependence of Eq. (|6|) allows for agreement with the 
COBE constraint with only a relatively small M variation. This is very promising from 
the perspective of relating inflation models to real scales of particle physics. To answer the 
questions of how well these ideas really work, and how constrained the parameters of the 
models really are, requires a detailed investigation of particular examples of these ideas. 



III. SUPERNATURAL INFLATION 

We define Flat Direction Hybrid Inflation (FDHI) models as those motivated by the 
properties of moduli fields or flat directions of the standard model. 

In the first model, we assume the existence of a superpotential which couples ip and 4> 
but which is suppressed by a large mass scale M'. For standard model flat directions, such 
higher dimension operators are to be expected, with M 1 equal to M p , M@, or some dynamical 
scale. In the case of moduli fields, it might be that this scale is of dynamical origin; one can 
readily determine how the answer changes with the form of the superpotential and the size 
of the mass scale. The example we take is 

w = £77- . 7 
2AP v ; 

We now need to specify the form of the supersymmetry breaking potential. We assume 
both if} and have mass of order the soft SUSY breaking scale of order 1 TeV (where we 
will need to test the consistency of this assumption). We assume that the potential for 
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the if) field gives a positive mass squared at the origin, while the field has negative mass 
squared at the origin. Furthermore, we assume that the cosmological constant is zero at the 
minimum of both if) and 0. The specific form of the potential we choose is 

r = «W(,/,/v) + ^ + **^, (8) 

where we have taken the scalar field to be real. When the parameters are motivated by 
supersymmetry breaking, we refer to our models by the name supernatural inflation. We 
will see that one very naturally obtains the correct magnitude of density perturbations, and 
sufficiently many e-foldings of inflation, using parameters and a potential which are well 
motivated in supersymmetric models. 

The if) mass is and the magnitude of the (imaginary) mass term (at the origin) is 
= M 2 1 f. During inflation, is confined near the origin. The field if> slowly rolls towards 

the origin and inflation ends about when if) — ip c — ^jlM'm^. It will turn out that either 
m^/m^ or M'/Mp is small, so that during inflation the term m^if) 2 is small relative to M 4 . 
The Hubble parameter during inflation is therefore approximately H = ^8ir/3M 2 /M p . 

We expect / is of order M P , or equivalently, m 2 ^ is of order m^. Although it looks like 
we took a very special form for the 0-potential, the use of the cosine is not essential. As can 
be seen from a Taylor expansion, only at the very late stages of inflation are terms other 
than the constant and mass term relevant. We could equally well have specified a potential 
which is truncated at fifth order in the fields, or which has different higher order terms. 
Although both if) and might be moduli or standard model flat direction fields, we assume 
their potentials are of very different form; the particular case we assume is illustrative of 
how a model could work. 

The constraint from density perturbations in the slow-roll regime is [||J| 
V 3/2 /[M^(dV/dip)] = 6 x 1(T 4 , where M p = Mp/y^n. This gives the constraint 

-^\lZl^ %m = 6-7 x 1(T 6 . (9) 

Here we have defined = m^/H and have measured time in e-foldings away from the 
time N = when if) — if) c (where inflation ends at positive N). It is clear that a lower 
M 1 makes the value of if) at the end of inflation lower, which in turn increases the density 
perturbations. The exponential in Eq. (Q) determines the scale dependence of the density 
perturbations, characterized by the scalar index. 

The scalar index a s is readily determined from the scale dependence of the density per- 
turbations to be This can be seen directly from the formula for density perturbations 
above. Alternatively, it is extracted from the general formula [H| 

(V'\ 2 V" 

n=l-2a s = l-3\ — \ +2—, (10) 

where dimensionful factors should be compensated by M p . Notice that the second term is 
negligible for all models for which the inflaton field value is much less than M p . This is 
readily seen from the fourth expression in Eq. @, which implies M P V /V ~ A<f>/(NM P ). 
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The third term is positive in our model, because the inflaton field rolls toward, rather than 
away, from the origin during the end of inflation. 

We see for this model that n is always greater than 1, and is very close to 1 for small 
/i^, which is the case for large M' . This differs from the usual prediction for new inflation 
or chaotic inflation models. The current upper bound on n is uncertain as is summarized 



in Ref. W% - These bounds, along with the validity of slow roll, prevent too large values 



of \x^. It turns out that in this model, the constraint for the correct magnitude of density 
fluctuations requires that when M' « M p , ~ 1/100, when M' ~ M G , /i^ ~ 1/10, and 
when M' « Mj, jj,^ ~ 1. These constraints are given in more detail in Ref. |12[| . 

Another distinctive feature of these models is that the ratio of the tensor to scalar 
contribution to the quadrupole 1Z = T/S ~ (V'/V) 2 ~ 0. Again this follows from the small 
value of the inflaton field if) near the end of inflation. 

As we have argued in the first section, models of inflation which have only a single field 
should have the inflaton field taking a value of order M p near the end of inflation if 50 
e-foldings are to be obtained without fine tuning. The combination of negligible 1Z and n 
never below 1 are distinctive features of these models which should help distinguish them 
from other possible inflationary models in the future. 

Another distinctive feature of the perturbation spectrum from this model will be a spike 
at small wavelengths. This spike can be calculated from a detailed study of the evolution of 



the if) and fields, as is performed in Ref. |12| . Here we give the qualitative features of the 
evolution. 

As we have emphasized, the if> field acts as a trigger to end inflation. While the field 
is confined to the origin, the evolution of the if) field is straightforward; it moves in towards 
the origin according to the slow roll equation of motion. 

While the if) field is big, the field has a large mass, so the field is strongly confined to 
the origin. As if) rolls in towards the origin, the <f> mass is decreasing. Once the <f> mass is 
sufficiently small, it random walks away from the origin due to de Sitter fluctuations. As 
the mass squared of the <f> field becomes negative, the field value becomes sufficiently large 
that the classical potential dominates over the de Sitter fluctuations as a driving force. The 
field then moves according to its classical equations of motion towards its true minimum. 

Because of the coupling of the if) and fields, the mass of the if> field increases as 
increases. Eventually, the time dependent if) mass is sufficiently large that the if) field begins 
to oscillate about its true minimum as a coherent state, and subsequently decay. Once the 
field oscillates about its true minimum, the exponential expansion of the universe ceases, and 
ultimately the field decays, permitting the universe to reheat to a temperature somewhat 
above the weak scale. 

The spike in the spectrum arises from fluctuations in the field during the stage in its 
evolution where de Sitter fluctuations drive its motion. We have calculated this spike fl2| 
and constrain our model so that inflation ends sufficiently rapidly that the spike is only 
relevant to perturbations on small (as yet unobservable) scales, less than 1 Mpc. This has 
given us a lower bound on fi^. We have also checked that the magnitude of the spike is 
consistent with known black hole constraints on these small scales. 

An alternative model can be constructed based on a renormalizable potential. Although 
this seems strange since we are considering flat direction fields, it is often the case, even in 
the MSSM, that fields cannot be simultaneously flat. So we take the potential to contain the 
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soft supersymmetry breaking terms as before but to contain a Yukawa coupling involving ip 
and 0. Specifically 



K=MW(0/V2/) + 2f! + A^. (11) 

This model has the essential features of the FDHI model of the previous section. The 
difference is the value of ip c which in this model is ip c = \/2m^/A. The density fluctuations 
give the constraint 

= = 1.6 x 1(T 4 . 12 

If A is of order unity, to satisfy Eq. (0) requires /x^ > 10 3 . However, if A ~ 10~ 4 — 10~ 5 , 
the model works perfectly with /j^ and both of order unity. We see there is virtually 



no fine tuning, so long as a small A exists. In Ref. |L2] we discuss examples, and construct 
an explicit GUT extension of the MSSM in which the small Yukawa coupling required by 
the particle physics provides the small coupling needed for the correct amplitude of density 
fluctuations. 



IV. CONCLUSIONS 

We have shown that with more than one field it is possible to construct models of inflation 
with no small parameters. Furthermore, the mass scales which seem to most naturally 
appear in these models are of order m 3 / 2 , about 1 TeV, and Mj, about 10 11 GeV, leading to 
a natural association with supersymmetric models. These models give rise to the correctly 
normalized density perturbations, even though the Hubble constant is quite low, of order 
10 3 — 10 4 GeV, because the value of the inflaton field at the end of inflation is much lower 
than the Planck scale. The key to producing more such models is a sensitive dependence of 
the (ft potential on the value of the ip field, so that the motion of the if) field can trigger the 
end of inflation while its value is small. 

An important feature which distinguishes these models from previous examples of hybrid 
inflation is that both the ip and <fi fields are assumed to have small mass. This is the 
more natural choice if both fields correspond to flat direction fields. Furthermore, because 
supersymmetry is broken during inflation, it is difficult to maintain the hierarchy between 
the two fields due to large radiative corrections. One would need further complications to 
maintain this hierarchy naturally. The small mass of <fi leads to the characteristic spike on 
small length scales. 

Our model is also consistent with the gravitino bound and Affleck Dine baryogenesis |12| . 



It seems that multifield models are probably the most natural models which can imple- 
ment inflation with weak scale Hubble constant, and that furthermore, these are probably 
the most natural inflation models in that they involve no new small parameters. The requi- 
site small parameters arise naturally from the ratio of mass scales. These models have the 
further advantages that they can be explicitly realized and one can calculate the relevant 
parameters for any particular implementation. They might even occur in simple extensions 
of the MSSM. 



S 



Perhaps the most important property of a model is its testability, and our proposed 
models have several characteristics that are in principle observable. The scalar index n 
which characterizes the scale dependence of density perturbations is always greater than 
unity. It is very close to unity for the model of Eq. (8) with M' at the Planck or GUT scale, 
but for M' at the intermediate scale or for the model of Eq. (11), it could be as large as 1.2 
for the parameters presented in our plots. In all cases tensor perturbations are negligible. 
An especially distinctive feature is a large spike in the density perturbation spectrum at 
present wavelengths of about 1 Mpc or less. 
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